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Figure 1: Standard Model contribution to K°-K° mixing (left) and a possible supersymmetric contribution 
(right). 



1. Introduction 

The Standard Model of particle physics has proven to be a remarkably successful theory; it 
describes the observed elementary particles and their interactions and explains the outcomes of 
high-energy collider experiments. Nevertheless, we know that the Standard Model is not a "theory 
of everything" because it does not account for neutrino masses, the amount of dark matter and 
dark energy in the universe, or the abundance of matter over antimatter. Nor does it explain the 
origin of particle masses and mixings. Specifically, in the quark flavor sector, we do not know why 
there are three generations of quarks, or what generates the hierarchy of quark masses and CKM 
mixing matrix elements. Thus we expect that the Standard Model is only a low-energy effective 
theory, and that we will discover new physics at higher energies that accounts for phenomena such 
as dark matter and explains the values of Standard Model parameters such as the quark masses and 
mixings. 

Most extensions of the Standard Model contain new CP- violating phases and new quark flavor- 
changing interactions. We therefore expect new physics effects in the quark flavor sector. The 
flavor sector is sensitive to physics at very high scales because new particles will typically appear 
in loop-level processes. For example, Fig. 1 shows a sample supersymmetric contribution to neutral 
kaon mixing. Thus we may see evidence for new physics in the flavor sector before we produce 
non-Standard Model particles directly at the LHC. 

Flavor factories such as those at BEPC, Cornell, Frascati, KEK, and SLAC have been pouring 
out data to measure hadron masses and mixings, branching fractions, and other observables. These 
results can be used to obtain CKM matrix elements, in some cases with percent-level accuracy. We 
can also compare the measured results with Standard Model predictions; significant inconsistencies 
would indicate the presence of new physics. Lattice QCD calculations are needed, however, to 
interpret many of the results of flavor physics experiments and to obtain the values of the CKM 
matrix elements. This is because in order to accurately describe weak interactions involving quarks, 
one must include effects of confining quarks into hadrons. Typically the nonperturbative QCD 
effects are absorbed into quantities such as decay constants, form factors, and bag-parameters that 
are computed in lattice QCD. 

In this review I present the status of lattice QCD calculations needed to determine CKM matrix 
elements and constrain the apex of the CKM unitarity triangle. Towards this aim, I highlight quan- 
tities either where significant computational progress has been made in the past few years, or where 
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there is an apparent tension with determinations from independent (non-lattice) methods. Although 
I intentionally omit most details about the lattice calculations, I provide references throughout this 
work for the interested reader. I also discuss some additional lattice QCD calculations that are of 
interest to flavor physics phenomenology, but for the most part have not been attempted. I focus on 
calculations of interest to rare B- and AT-decays. I emphasize, however, that this list is by no means 
exhaustive, and that more lattice QCD calculations of weak matrix elements will become possible 
as computational resources increase and lattice QCD methods improve. 

This paper is organized as follows. First, in Sec. 2, 1 briefly review the Standard Model CKM 
framework. Next, in Sec. 3, I summarize the current status of the CKM matrix, focusing on de- 
terminations of CKM matrix elements and constraints on the CKM unitarity triangle that rely on 
inputs from lattice QCD. In the first subsection, Sec. 3.1, 1 present determinations of the elements 
of the first row of the CKM matrix; the values are currently consistent with the Standard Model 
expectation of unitarity. In the second subsection, Sec. 3.2, I present current averages for lattice 
QCD inputs needed to constrain the CKM unitarity triangle. I also show the potential impact of im- 
proved lattice QCD calculations of these quantities on the global CKM unitarity triangle fit. Next, 
in Sec. 4, 1 discuss several hints of new physics in the flavor sector that rely on lattice QCD calcu- 
lations of weak matrix elements. For example, given current theoretical and experimental inputs to 
the CKM unitarity triangle analysis, there is a ^ 2-3 a tension between observations and Standard 
Model expectations. In Sec. 5, I discuss prospects for lattice QCD calculations to search for evi- 
dence of new physics in rare B- and /^-decays. Finally, I conclude by discussing future prospects 
for lattice QCD to aid in new physics searches in the quark flavor sector in Sec. 6. 



2. The Standard Model CKM framework 



In this section I briefly review the Standard Model CKM framework. For a more thorough 
discussion, see, for example, Ref. [1]. 

The 3x3 Cabibbo-Kobayashi-Maskawa (CKM) matrix parameterizes the mixing between 
quark flavors under weak interactions [2, 3]: 
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where the numerical values are from Ref. [4]. The matrix elements are empirically largest along 
the diagonal, so mixing is most probable within the same generation. Because the CKM matrix is 
unitary in the Standard Model, one can express elements of the CKM matrix as an expansion in 
powers of the small parameter X = \V US \ ~ 0.22. The resulting Wolfenstein parameterization [5], 
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Figure 2: "Gold-plated" processes on the lattice that can be used to obtain each CKM matrix element. 

n — 

Neutral K v -K mixing is also a gold-plated process, and gives a constraint on the phase of the CKM matrix 



(p/n) 




(0,0) (1,0) 
Figure 3: The CKM unitarity triangle. Figure from Ref. [4]. 

makes the hierarchy of sizes among the matrix elements explicit in terms of powers A. 

Because the CKM matrix elements are fundamental parameters of the Standard Model, it 
is important to know their values as precisely as possible. "Gold-plated" lattice processes allow 
the determination of most CKM matrix elements, with the exception of V t b- These are simple 
processes with only one hadron in the initial state and at most one hadron in the final state, where 
the hadrons are stable (or at least narrow and far from threshold). They are therefore the easiest 
to calculate with standard lattice QCD methods, and can now be obtained reliably from "2+1" 
flavor simulations that account for the dynamical w-, d-, and s-quarks with all sources of systematic 
error under control [6]. Figure 2 lists the gold-plated lattice quantities that can be used to obtain 
each CKM matrix element. For a review of lattice calculations of most of these quantities, see 
the plenary talks in these proceedings by C. Aubin on heavy quark physics [7], V. Lubicz on kaon 
physics [8], and E. Scholz on light pseudoscalar meson masses and decay constants [9]. 

The quark mixing matrix is unitary, leading to normalization and orthogonality relations among 
the matrix elements. In the case of three generations of Standard Model quarks, the most important 
relationship is 

V ud V: b + V cd V c *„ + V td V t *„ = . (2.4) 

Equation (2.4) can be expressed as a triangle in the complex p-T] plane known as the CKM unitarity 
triangle; this is shown in Fig. 3. Note that the CKM unitarity triangle is rescaled by IV^V^I so 
its base has unit length. In practice, terms of and higher must be included in order to achieve 
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sufficient accuracy for phenomenology. Thus it is standard to work with the parameters p and rj, 
which are defined through all-orders in X as the apex of the triangle in Fig. 3: 



p + irj = — 




(2.5) 



Constraints on (p , Y\) can translated in a straightforward manner into constraints on original Wolfen- 
stein parameters (p , T] ) . 

Deviations from the relationship in Eq. (2.4) would indicate the presence of physics beyond the 
Standard Model. It is likely that new physics will have additional quark flavor changing interactions 
and CP- violating phases. Once theory and experiment are sufficiently precise, these will manifest 
themselves as inconsistent determinations of the apex of the unitarity triangle (p, T]). One of the 
key goals in flavor physics is therefore to determine (p , ff) as precisely as possible from as many 
independent processes as possible in order to search for new physics. 

Several important constraints on the CKM unitarity triangle require both experimental mea- 
surements and lattice QCD calculations of nonperturbative hadronic weak matrix elements. Lattice 
calculations of neutral kaon mixing are needed to determine the mixing parameter B K and constrain 
the phase of the CKM matrix. Computations of neutral 5-meson decays and mixing are needed to 
obtain the decay constants fs and mixing parameters B Bq (q = d,s)\ these are used to determine 
\V t d\ and \V ts \. Calculations of the semileptonic 5-meson decay B — ► ntv are needed to compute 
the form factor f+(q 2 ), which allows a determination of \V u t\, while those of the decays B — > D£v 
and B ^D*£v are needed to compute the form factors F ( 1 ) , G( 1 ) , which allow one to obtain \V c b\> 

3. Status of the CKM matrix 

In this section I focus on tests of the Standard Model CKM framework, emphasizing the role 
played by lattice QCD inputs. In Sec. 3.1 1 review the status of first-row unitarity, while in Sec. 3.2 
I review that of the global CKM unitarity triangle fit. 

3.1 First-row unitarity 

In the Standard Model, elements of the first row of the CKM matrix must obey the following 
relation: 

\Vud\ 2 + \V U s\ 2 + \Vub\ 2 = l. (3.1) 

Because \V u b\ ~ <^(10 -3 ), Eq. (3.1) is essentially a constraint on the relationship between \V u d\ and 
\V US \. A significant deviation in this relationship would be evidence for new physics. 
The CKM matrix elements \V u d\ and \V US \ are known to sub-percent accuracy [10]: 

\V ud \ = 0.97418 ±0.00026, (3.2) 
\V US \ = 0.2246 ±0.0012, (3.3) 

where \V u d\ is obtained from nuclear j8 -decays [11] and \V US \ relies on a Nf = 2+ 1 flavor lattice 
QCD calculation of the K —> ntv form factor by the RBC and UKQCD collaborations [12]. The 
ratio iKsl/lKdl is also known to this accuracy [10]: 

\V us \/\Vud\= 0.2321 ±0015, (3.4) 
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Figure 4: Combined fit to \V u d\, \V US \, and IV^I/lK^I both with and without the constraint from first-row 
unitarity [10]. The fit with first-row unitarity imposed has a # 2 /d.o.f. = 0.65 and C.L. = 42%. 

which uses a Nf = 2 + 1 flavor lattice QCD determination of fx/fn by the HPQCD collabora- 
tion [13]. Equation (3.1) therefore allows for a precision test of the Standard Model and probe of 
new physics. 

The Flavianet collaboration has performed a combined fit to \V u d\, \V US \, and IV^I/lKdl using 
the values in Eqs. (3.3) and (3.4), both with and without imposing the constraint from first-row 
unitarity [10]. This is shown in Fig. 4. Without assuming unitarity, the result of the combined fit is 

1^1=0.97417(26), (3.5) 
1^1=0.2253(9), (3.6) 

which is consistent with first-row unitarity at subpercent level: 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.9998(6) . (3.7) 

In addition, the fit including the unitarity constraint has a confidence level of 67%. Thus there is 
no indication of new physics in the first row of the CKM unitarity triangle given the current levels 
of experimental and theoretical precision. 

3.2 The CKM unitarity triangle 

In this section I summarize the current status of the global fit of the CKM unitarity triangle. For 
this discussion, I use the recent averages of the necessary lattice QCD inputs provided in Ref. [14]; 
these are given in Table 1. These results were obtained by averaging all available Nf = 2+1 
flavor results documented in proceedings in publications that contain complete error budgets, and 
therefore do not include any new results as of Lattice 2009. See the plenary talks of C. Aubin [7], 
V. Lubicz [8], and E. Scholz [9] in these proceedings for discussions of more recent preliminary 
and Nf = 2 results. For the computation of the lattice QCD averages in Ref. [14], asymmetric 
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Quantity Error 

B K = 0.725 ±0.026 -4% 

/^a/^ = 275±13 -5% 

§ = 1.243 ±0.028 -2% 

iK^lexci. = (3.42±0.37) x 10- 3 - 11% 

iK^lexci. = (38.6 ± 1.2) x 10~ 3 - 3% 

Kg =0.92±0.01 - 1% 

/* = 155.8 ±1.7 MeV - 1% 



Table 1: Current status of lattice inputs to the global fit of the CKM unitarity triangle [14]. These were 
obtained by averaging all available Nf = 2 + 1 results documented in proceedings and publications that 
contain complete error budgets, and account for correlations between different calculations in a conservative 
manner. 



HFAG Winter '09 (no B->sy) 
exclusive B->Dlv 

exclusive B->D lv 

37 38 39 40 41 42 43 
IV J X 10" 3 

CD 

Figure 5: Tension between inclusive (blue square) and exclusive (red circles) determinations of \V c b\. The 
inner and outer errors bars on the inclusive value denote the errors from the moment fit and from other 
systematics, respectively. The inner and outer errors on the exclusive values are experimental and theoretical, 
respectively. The inclusive determination [15] disagrees with the exclusive determination from B — > D£v by 
1.2<7 and from B D*tv by 2.3d. 

errors are symmetrized and it is assumed that errors are normally distributed. Furthermore, when 
there are correlations between a source of error in different lattice calculations, it is conservatively 
assumed that the correlation is 100%. 

3.2.1 iKalandlV^I 

Tensions exist between the inclusive and exclusive determinations of both \V c b\ and \V u b\. The 
inclusive values are independent of lattice QCD, while the exclusive results rely on lattice compu- 
tations of 5-meson semileptonic form factors. 

Figure 5 compares the inclusive determination of \V c t\ [15] with exclusive determinations 
based on lattice QCD calculations of the B — > D£v [16] and B — > D*£v [17] form factors at zero 
recoil by the Fermilab Lattice and MILC collaborations. The tension between inclusive and exclu- 
sive values is ~ 1.2-2.3(7. This discrepancy is significant because \V c b\ is used to normalize the 
base of the unitarity triangle. In fact, the error in \V c b\ is now the limiting uncertainty in some of 
the unitarity triangle constraints, as shown in the next subsection. Therefore improving the errors 
in exclusive \V c b\ should be a high priority for lattice QCD calculations. 
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Figure 6: Tension between inclusive (blue squares) and exclusive (red circle) determinations of \V u b\- The 
inner and outer errors on the inclusive values are experimental and theoretical, respectively. The error bar 
on the exclusive value denotes the total statistical plus systematic error added in quadrature. The exclusive 
determination disagrees with the inclusive determinations [15] by 1-2 a depending upon the theoretical 
framework used to obtain the inclusive value. 



Figure 6 compares two representative inclusive determinations of \V u b\ [15] with the exclu- 
sive determination based on lattice QCD calculations of the B — > n£v form factor by the HPQCD 
collaboration [18] and by the Fermilab Lattice and MILC collaborations [19]. The l-2<7 tension 
between inclusive and exclusive determinations of \V u b\ is less worrisome, however, than in the 
case of \V c b\ because the inclusive value of \V u b\ varies significantly depending upon theoretical 
framework, and is highly sensitive to the input value of the Z?-quark mass [20]. Nevertheless, if the 
discrepancy in \V u b\ holds up as the size of the theoretical errors are reduced, it could be a sign of 
non-Standard Model (V + A) currents. A right-handed b —> u current would affect the extraction 
of \V u b\ from B — > n£v semileptonic decay differently than the extraction of \V u b\ from B — > X u £v 
inclusive semileptonic decay because B — > n£v decay proceeds through only the vector part of the 
left-handed current, whereas B — » X u £v inclusive decay is the sum of several exclusive channels that 
proceed via different linear combinations of the vector and axial- vector currents. The right-handed 
current hypothesis could be tested using \V u b\ determined from exclusive B — > p£v semileptonic 
decay since B — > p£v proceeds through the entire left-handed current. One could use the values of 
\V u b\ from B — > 7t£v and B — > X u £v to make a prediction for the size of the non-Standard Model 
(V +A) current, and make a prediction for the value of \V u b\ from B — ► p£v given this hypothesis. 
This motivates the need for a reliable lattice QCD calculation of the B — ► p£v form factors. 

3.2.2 Neutral kaon mixing 

The amount of direct CP-violation in the neutral kaon system, £ K , which has been measured 
experimentally to sub-percent precision, constrains the apex of the CKM unitarity triangle: 



Until recently, the uncertainty in the unitarity triangle constraint from e K was primarily due to the 
uncertainty in lattice QCD calculations of the hadronic matrix element Bk- This is no longer true, 
however, thanks to the Nf = 2 + 1 calculations by the RBC and UKQCD collaborations [21] and by 
Aubin, Laiho, and Van de Water [22], which constrain Bk to ~ (4 — 6)% accuracy. Furthermore, the 
errors in B K will continue to decrease as the RBC and UKQCD collaborations update their analysis 



\e K \ = K £ C e £M 2 7H-mSo(*c)(l-^^ 



(3.8) 
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using two lattice spacings [23] and as independent determinations using twisted-mass fermions [24] 
and improved staggered fermions [25-28] become available. 

The dominant error in the £k band is now from the Wolfenstein parameter A, which is pro- 
portional to the CKM matrix element \V c b\. Although \V c b\ is known to ~ 2% accuracy, it enters 
Eq. (3.8) as the fourth power and therefore contributes an ~ 8% error to the £k band. Figure 7 
illustrates this point by showing the relative sizes of the contributions from B K and \V c b\ to the £k 
band. It is easy to see that further improvements in B K will not tighten the constraint from e K on 
the apex of the CKM unitarity triangle without an improvement in \V c b\. 

Given the recent reductions in the errors on Bk and \V c b\, it is now important to include 5 — 10% 
corrections to Eq. (3.8) that have been omitted from most previous unitarity triangle analyses. In 
the Standard Model, direct CP- violation is parameterized as: 



where the first term in parentheses is the short-distance contribution which is proportional to B K 
and the second term, Po, is the long-distance contribution which is related to the ratio of K — > kk 
decay amplitudes in the AI = 1/2 channel. Because the phase (f) £ is close to 45° and the size of 
Po is small compared to the short-distance contribution, it is common to set (f) £ = 45° and Po = 0; 
this leads to the standard expression given in Eq. (3.8). Buras and Guadagnoli recently pointed out, 
however, that the corrections due to (f) £ ^ 45° and Po ^ 0, are no longer negligible given the sizes 
of other errors [29]. 

The effects of (j) £ ^45° and Po ^ can be parameterized by an overall multiplicative fac- 
tor [30]: 



where K £ is given by 



(3-10) 

AM K 

K £ = Vising ( 1 - -W/ £ ,) + h ^^) • (3.11) 
V w V2\e K \ J 
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Figure 8: Global fit of the CKM unitarity triangle [14] . The current fit is consistent with the Standard Model 
at the 23% level. The constraints from |K& 1/1^61 > AM^/AMj, and AMj are all limited by theoretical 
uncertainties from lattice QCD. 

Fortunately, assuming the Standard Model, the majority of inputs to K £ are well-known from exper- 
iment. The remaining unknown, ImA2, can be obtained from lattice QCD calculations of K — > kk 
matrix elements in the AI = 1/2 channel. The only Nf = 2+1 flavor determination of this quantity 
is by the RBC and UKQCD collaborations [31], and has rather large systematic errors associated 
with the use of leading-order chiral perturbation theory. Nevertheless, the contribution to Eq. (3.11) 
from ImA2 is small and leads to only a ~ 1% uncertainty in K e [14]: 



This result agrees with the estimate of Buras and Guadagnoli [29] and lowers the SM prediction 
for £k by 8%. The correction factor K e is included in the global unitarity triangle fits presented in 
the following section [14], and has also recently been included by the UTfit collaboration [32]. It 
has not yet been implemented by the CKMfitter group [33]. 

3.2.3 Global fit of the CKM unitarity triangle 

Figure 8 shows the current status of the global CKM unitarity triangle fit using the lattice QCD 
inputs presented in Table 1 [14] . Although the average of inclusive and exclusive determinations 
of \V c b\ is used, the error in the average is inflated in order to account for the inconsistency be- 
tween the two values following the prescription of the Particle Data Group [4]. Only the exclusive 
determination of \V u b\ is used, however, because the inclusive determination varies so much de- 
pending on the theoretical framework. The confidence level of the global fit is 17%; thus current 
observations are consistent with Standard Model expectations given the present level of theoretical 
precision. 

Currently the constraints from £ K , Am^/Amj, and \V u b/V c b\ are limited by uncertainties in the 
lattice QCD calculations of IV^Iexci.* an d |Kfr|exci.> respectively. In order to show the poten- 
tial impact of future lattice calculations, it is therefore an interesting exercise to repeat the global 
CKM unitarity triangle fit after reducing the lattice uncertainties in Bk, iKfclexci.* an d |Kfr|exci. 
to 1% (with central values fixed). The resulting fit is shown in Fig. 9. In this case, only the exclu- 
sive determination of \V c b\ is used because combining it with the inclusive determination becomes 



Ke = 0.92 ±0.01. 



(3.12) 
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Figure 9: Potential impact of future lattice determinations on the global unitarity triangle fit. If the theoret- 
ical errors in all of the lattice QCD inputs are reduced to 1 % with the central values fixed, the fit would no 
longer be consistent with Standard Model expectations. Figure courtesy of E. Lunghi. 

problematic when the lattice errors are so small. The confidence level of this hypothetical fit is 
0.0004%, indicating that this scenario would no longer be compatible with the Standard Model. 
Thus improved lattice QCD calculations of hadronic weak matrix elements could allow the obser- 
vation of new physics in the quark flavor sector with a high significance. 

4. Hints of new physics in the flavor sector 

Although most observations in the flavor sector are consistent with Standard Model expecta- 
tions, there are currently several 2-3(7 tensions that may indicate the presence of new physics. In 
the following subsections I present those hints of new physics that rely on lattice QCD calculations 
of weak matrix elements. It is worthwhile, however, briefly mentioning first two others that do not 
require lattice inputs. In the case of all tensions, we must, of course, wait and see whether their 
significance increases or decreases with improved experimental and theoretical precision. 

The -mixing phase in the Standard Model is given by j8, SM = arg(-V ts V* b /V cs V* b ) « 0.02. 
This prediction disagrees, however, with the measured world average based on flavor-tagged anal- 
yses of B s — > J/y/cj) decays by 2.2(7 [34]. If this discrepancy remains as the experimental errors 
are reduced, it would indicate the presence of new sources of CP-violation beyond the phase of 
the CKM matrix. Such new phases would lead to correlated effects between AS = 2 processes and 
b — > s decays, thereby making improved measurements of CP- violation in b — > s penguin decays 
particularly important [35]. Thus is is interesting that the "effective" value of the angle j8 in the 
CKM unitarity triangle obtained from b — » qqs penguin decays is lower than the value of j8 ob- 
tained from tree-level b — > ccs decays, which are expected to be less sensitive to new physics. For 
example, sin(2j8) e ff determined from the penguin decay B — » <pK° is ~ 1.3(7 lower than the average 
sin(2j8) from tree-level decays [36]. Although this discrepancy is not statistically significant, the 
value of sin(2j8) e ff obtained from various penguin decay modes is systematically lower than the 
average sin(2j8) from tree-level decays in almost all decay channels. Addressing these puzzles will 
require improved determinations of the fi^ -mixing phase and of sin(2j8) e ff from b — ► qqs penguin 
decays at LHCb and the super-2? factory. 
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Figure 10: Tension in the CKM unitarity triangle [14]. The confidence levels shown in the figures corre- 
spond to the unitarity triangle fit including all three constraints, while the solid, dashed and dotted contours 
are obtained by omitting S^k and AMbJ AMs d , respectively. The left and right panels use the exclusive 
and inclusive determinations of \V c b\- 



4.1 Tension in the CKM unitarity triangle 

As noted in the introduction, there is a 2-3(7 tension in the CKM unitarity triangle [37, 32, 14]. 
If the constraints on (p , f[ ) from a, y, and | V u b | are omitted, it is easy to see that the tension is really 
between the three most precise constraints, which are those from £k, Am^/Amj, and sin(2j8) [37]; 
this is illustrated in Fig. 10. Furthermore, the degree of tension is extremely sensitive to the value 
of \V c b\> The confidence level of a fit to £k, Am^/Amj, and sin(2j8) using the exclusive determi- 
nation of \V c b\ is only 0.2% (left plot of Fig. 10), whereas the confidence level using the inclusive 
determination is almost 50 times greater at 8.9% (right plot of Fig. 10). 

One way to quantify the significance of the tension is to leave one input as a free parameter 
in the unitarity triangle fit and make a prediction for it based on the remaining constraints. For 
example, the value of \V c b\ preferred by £k, sin(2/3) and AM Bs /AM Bd is [14]: 

|^| fit = (43.0±0.9)xl0- 3 . (4.1) 

This deviates by 3.0(7 and 1.3<7 from the exclusive and inclusive determinations of \V c b\, respec- 
tively. Because of the sensitivity of the unitarity triangle analysis to the value of IV^I (and the 
discrepancy between inclusive and exclusive determinations), obtaining \V c b\ should be a high pri- 
ority for lattice QCD calculations. 

Alternatively, the value of B K preferred by the other inputs depends upon the value of | V c b | [14]: 

fl.09db0.12 |V cfe | excl 
(B K U= < 0.903 ±0.086 |V cfo | incl (4.2) 
[ 0.98 ±0.10 |V c ,| excl+incl 

These deviate from the lattice QCD average for Bk at the 2.9(7, 2.0(7 and 2.4(7 levels, respectively. 
The observed tension between the £k band and the other unitarity triangle constraints relies on the 
inclusion of the correction factor K £ , which raises location of the £k band in the ~p-ff plane by 8%. 

4.2 The "f Ds puzzle" 

The HPQCD collaboration's published lattice QCD calculation of the leptonic decay constant 
fu s [13] disagrees with the CLEO experimental measurement [38] by 2.4(7, despite the fact that 
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Figure 11: Comparison of experimental measurements and lattice QCD calculations of fo and fjy s . The 
majority of the results shown agree, but the HPQCD calculation of fn s disagrees with the CLEO measure- 
ment by 2.4cr. This discrepancy has yet to be confirmed, however, by an independent lattice calculation with 
comparably small theoretical errors. Furthermore, the HPQCD determination of fo s will likely increase by 
~ 1.5cr when the recently updated value of the scale r\ [45] is taken into account. 



HPQCD's determinations of f K , /k, /d all agree with experiment. Figure 11 shows a compar- 
ison of recent experimental and theoretical determinations of fr> and fo s . This disagreement is 
somewhat challenging to accommodate in new physics models because Standard Model D s lep- 
tonic decay occurs at tree-level. Dobrescu and Kronfeld have pointed out, however, that models 
with a charged Higgs or leptoquark can work [39], and that these could also lead to a signal in 
D — » K£v semileptonic decay [40]. The only other Nf = 2+1 flavor calculation of fo s by the 
Fermilab Lattice and MILC collaborations, which was updated in these proceedings [41], agrees 
with the CLEO measurement, but has an error that is approximately three times larger than that of 
HPQCD. The ETM collaboration recently published an Nf = 2 determination of fo s using twisted- 
mass fermions that agrees with the CLEO measurement [42], but the result also has larger errors 
than those of HPQCD. Alternative lattice calculations of and fp s with comparable precision to 
that of HPQCD are therefore necessary to either confirm or refute the result. On the experimental 
side, BES-III should measure the D- and ZVmeson leptonic decay constants with ~1% precision 
after 4 years of running at their nominal luminosity [43], and will help shed light on this 2.4(7 
tension. 

It is worth noting that the HPQCD collaboration recently obtained a new value for the absolute 
scale r\ which they use to convert lattice quantities into physical units. Their published result 
for f Ds relies upon the older value of r\ = 0.321(5) fm from the T 1S-2S mass-splitting [44]. 
Their new r\ = 0.3133(23)(3) fm is ~ 1.5(7 lower, and comes from averaging the value of r\ 
obtained from three different methods [45]. Because the uncertainty in f Ds is currently dominated 
by the uncertainty in n, use of the new scale determination will likely increase the value of fo s 
by approximately 1.5<7, or 5 MeV, and thereby reduce the tension with the CLEO measurement to 
below 2.0(7. Thus the "fo s puzzle" may not, in fact, be a puzzle for much longer. 
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Figure 12: Tension between the Standard Model prediction for B — > TV leptonic decay and current experi- 
mental observations. The left plot compares lattice QCD calculations of fs with experimental determinations 
of fs (using \Vub\ from lattice QCD) [46]; the discrepancy is 1.9<7. The right plot compares the Standard 
Model prediction for the B — ► TV branching fraction obtained from the global unitarity triangle fit with the 
experimentally-measured world average [50]; the discrepancy is 2.4o\ 



4.3 B — > TV leptonic decay 

There is a 1.9a difference between lattice calculations of fs and experimental determinations 
of fs (using \Vub\ exclusive from lattice QCD), as shown in the left-hand plot in Fig. 12 [46]. 
This discrepancy is unrelated to the fr> s puzzle discussed in the previous subsection. The HPQCD 
collaboration's calculation of fs [47] uses NRQCD Z?-quarks, and is independent of their determi- 
nation of fo and fo s using HISQ charm quarks. Furthermore, this difference is also observed by 
the Fermilab Lattice and MILC collaborations, whose determination of fs has a smaller error than 
that of HPQCD [48]. It should be noted both the HPQCD and Fermilab/MILC determinations of 
/b will increase slightly when they are recomputed using updated values of the scale r\ [45, 49]. 
This should not, however, affect the tension significantly because the uncertainty in fs from r\ is 
much less than the total error. 

There is also a 2.4-2.5(7 discrepancy between the Standard Model prediction for the B — > TV 
branching fraction from the global unitarity triangle fit and the experimentally-measured world 
average [50, 51], shown in the right-hand plot in Fig. 12. This comparison relies on several lat- 
tice QCD inputs to the various unitarity triangle constraints. It can be made independent of the 
theoretical uncertainty in f B , however, by using the ratio of £%{B — > Tv)/Amj. The CKMfitter 
group uses this ratio to make a prediction for the 2^-meson mixing parameter B Bd . They obtain 
Bs d = 0.5ltoio [52], which differs considerably from the only Nf = 2+1 flavor lattice result 
B Bd = 1.26±0.11 [53]. 

The discrepancy inB — > TV leptonic decay could be due to the presence of a charged Higgs, but 
this is not a particularly natural explanation. Although an enhancement over the Standard Model 
branching fraction could occur for very light values of the charged Higgs mass, most MSSM sce- 
narios would lead to a suppression of the branching fraction [54, 1]. Reducing the experimental 
errors in 3&{B —> tv) will be difficult at LHCb, but prospects are good for a super-2? factory [55]. 
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Figure 13: Standard Model contribution to B — > Kit semileptonic decay. This process is suppressed in the 
Standard Model because it proceeds through flavor-changing neutral currents. 

The theoretical errors in fs will improve more quickly. The Fermilab Lattice and MILC collabo- 
rations are reducing their errors by adding statistics, while the HPQCD collaboration is exploring 
the use of HISQ light quarks in combination with NRQCD Z?-quarks [56]. Ongoing work by the 
RBC and UKQCD collaborations using domain-wall light quarks and static ^-quarks [57] and by 
the ETM collaboration using twisted-mass light quarks and static ^-quarks [58] was also presented 
at this conference. 

5. Lattice flavor physics beyond the CKM matrix 

Because current observations are consistent with the Standard Model CKM framework with 
only a few 2-3 a deviations that may or may not turn out to be significant when experimental and 
theoretical errors improve, it is clear that we should look for new physics in other realms of the 
quark flavor sector beyond the CKM matrix. One set of processes that are particularly sensitive to 
physics beyond the Standard Model are "rare" decays, which are rare because they are suppressed 
at lowest-order in the the Standard Model. This can be because the decay is loop-suppressed (as in 
flavor-changing neutral currents), the decay is helicity-suppressed by powers of small masses, or 
the decay is "Cabibbo"-suppressed by powers of small CKM matrix elements. 

The potential to search for new physics in rare decays is far from exhausted. The Super-5 
factory, "Project X" at Fermilab, and other future experiments will improve measurements of (or 
observe for the first time) many rare B- and ^-decays. In this section I therefore discuss ways in 
which the lattice community can aid in the search for new physics in rare decays. 

5.1 Rare b — ► d and b — » s transitions 

Semileptonic decays involving b —> d and b —> s transitions are rare because the lowest-order 
Standard Model contributions are 1-loop and hence suppressed. For example, Fig. 13 shows a 
sample Standard Model contribution to B — > Kit decay. Because new particles will enter b — » d 
and b — > s electroweak penguin loops in many models, new physics contributions to these processes 
could be of the same size or larger than those from the Standard Model. Thus b — > d and b — > s 
transitions are potentially stronger probes of new physics than b — > u transitions. Figure 14 shows 
the differential branching fractions for B —> K*£ + £~ and B — » K£ + £~ versus q 2 measured by the 
Belle collaboration [59]. Although the errors are still rather large, both LHCb and Super-2? will 
soon improve the measurement of these and other rare 5-decays. 
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q 2 (GeV 2 /c 2 ) 

Figure 14: Differential branching fractions for B — » K*£+£~ (upper plot) and B — > K£ + £~ (lower plot) 
versus q 2 [59]. The two shaded regions are veto windows to reject J/\j/ and yf(2S) events. The solid curves 
show the Standard Model predictions. The widths of the allowed ranges are dominated by the theoretical 
uncertainties in the hadronic form factors [61, 62]. 

The Standard Model predictions for the B — > KS*>tt branching fractions (denoted by the red 
solid curves in Fig. 14) are limited by the uncertainty in the hadronic form factors. Estimates from 
light-cone sum rules typically lead to ^30% hadronic uncertainties in the branching ratios [60]. 
For example, in the case of B — > K££, the Standard Model prediction is [62] 

3B{B^Ktr) = (0.35 ±0.11 ±0.04 ±0.02 ±0.0005) x 10~ 6 , (5.1) 

where the first error is the contribution from the hadronic form factors. Because the branching 
fraction is proportional to the square of the hadronic form factors, a reduction in the errors on 
f+{q 2 ) and fr(q 2 ) to below 10% would reduce the total theoretical error significantly. Lattice 
QCD therefore has a good opportunity to improve the theory of rare 5-decays in the Standard 
Model. Once the hadronic matrix elements are computed precisely, they can also be combined 
with the perturbatively-calculated Wilson coefficients to make predictions for beyond-the-Standard 
Model scenarios and help distinguish between new physics models. 

Currently there are no published Nf = 2+1 flavor lattice QCD calculations of semileptonic 
form factors for rare ^-decays, although work is in progress using staggered light quarks and mov- 
ing NRQCD heavy quarks [63, 64]. As Fig. 14 shows, in the case of B^K<-*ke the full a 2 range is 
experimentally accessible (with the exception of the // \jf and \jf(2S) peaks). Thus the fact that stan- 
dard lattice QCD methods can only obtain semileptonic form factors at high q 2 is less problematic 
than in the case of B —> n£v decay. Furthermore, model-independent formulas based on analyticity 
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and crossing symmetry such as the "z-parameterization" can be used to extrapolate lattice form fac- 
tor results to low-g 2 in a clean way [65-67]. This method has already been employed in the case 
of B — > ntv by the Fermilab Lattice and MILC collaborations [19]. A challenge for lattice calcu- 
lations of rare 5-decays is the fact that many interesting processes such as B — > K*U and B — ► K*y 
involve vector mesons. In these cases chiral perturbation theory cannot be used to guide the chiral 
extrapolation. Fortunately, this situation will improve as the simulated quark masses continue to 
become lighter. 

5.2 K — > 7rvv decays 

The semileptonic decay K — > nv v proceeds through flavor-changing neutral currents, as shown 
in Fig. 15, so its lowest-order contributions are 1-loop and hence small in the Standard Model. In 
fact, K — > 7TVV decays are so rare that, although the process has been observed in the K + chan- 
nel [68], only a limit has been set in the more challenging K° channel [69]: 

^(K + -> 7T + vv) exp . = 1.73+1;^ x KT 10 , (5.2) 
3B{jQ -> 7r°vv) exp . < 670 x 10" 10 (90%CL) . (5.3) 

Both of these results will improve, however, in the near future. Experiment NA62 at the CERN SPS 
will measure ^(100) events in the K + channel, while E391a at KEK (to become E14 at J-PARC) 
will collect the first K° events. Thus it is important to consider how lattice QCD calculations can 
aid in the search for new physics in rare /^-decays. 

The Standard Model branching fractions for K — > 7TVV are theoretically under much better 
control than in the case of rare 5-decays, and are known to 10-15% accuracy [70, 71]: 

SS{lQ -> 7T°vv)sm = (0.249 ±0.039) x lO" 10 , (5.4) 
38{K + -> tt + vv)sm = (0.85 ±0.07) x 10" 10 . (5.5) 

Furthermore, the uncertainties in the hadronic form factors lead to only a tiny percentage of the total 
error. Figure 16 shows the theoretical error budgets for the Standard Model branching fractions in 
the K° and K + channels; the form factor Kl contributes 4% of the error in — * 7T°vv), while 
the form factor K + contributes 2% of the error in £${K + — > 7T + vv). This is because the form 
factors can be estimated quite precisely using experimental K —> 7t£v data combined with NLO 
and partially-NNLO ^PT to account for isospin-breaking effects [70]. 
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Despite the fact that the K — > nvv hadronic form factors are already well-known, lattice QCD 
can still help reduce the theoretical uncertainties in the Standard Model branching fractions in 
several ways. The dominant uncertainties in the K° and K + branching fractions are parametric 
errors from the CKM matrix elements. Because the branching fractions are proportional to the 
Wolfenstein parameter A 4 , lattice QCD can make a significant impact by reducing the uncertainty 
\V c b\. The subleading errors in 38 (K® — > 7T°vv), however, are from the uncertainty in the top 
quark mass and unknown higher-order corrections to the Inami-Lim function X(x t ), and must be 
addressed with better measurements of m t and higher-order perturbative calculations of X(x t ). The 
second-largest error in 38 (K + — > 7T + vv) is due to the long-distance contribution from up and charm 
quarks. Isidori, Martinelli, and Turchetti have proposed a method for computing 8P C U using lattice 
QCD that can reduce the contribution to the error on the rate to the 1-2% level [73]. The error 
budget for 38 (K + — > 7T + vv) presented in Fig. 16 uses the charm quark mass recently obtained 
from lattice QCD calculations of current-current correlators by the HPQCD collaboration [74]; 
thus the 6% error can be reduced even further with improved lattice determinations of m c . The 
remaining errors in 38 {K + — » 7T + vv), however, can only be addressed with better measurements 
of m u more precise determinations of a s , and higher-order perturbative calculations of X(x t ). 

Once the K —> Ttvv branching fractions have been measured more precisely, they can be used 
to determine the apex of the CKM unitarity triangle: 



?(K + -> 7T+VV) = K + A A X(x t f 



3§{K° L -> 7T°vv) = K L A 4 Tj 2 X(x t ) 2 . (5.7) 

These constraints are shown in Fig. 17. The location of the apex of the CKM unitarity triangle as 
determined strictly from kaons can then be compared with one from clean 5-physics observables 
like sin(2j8) and neutral 5-mixing. This will provide a highly non-trivial test of the CKM picture. 
Even if the K- and 5-physics predictions agree, rare kaon and 5-decays still probe new physics 
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operators at effective scales up to several TeV or higher. Thus they will play an important role in 
discriminating between Standard Model extensions. 

6. Summary and outlook 

Lattice QCD and flavor physics are entering a mature era in which experimental measurements 
and theoretical calculations are sufficiently precise to constrain the presence of new physics in 
the flavor sector. In particular, lattice QCD can now be used reliably compute hadronic matrix 
elements that encode the nonperturbative QCD contributions to weak processes. For many of 
the weak matrix elements of interest, there are now or soon will be several independent Nf = 
2 or Nf = 2+ 1 flavor lattice calculations, thereby lending credibility to the results [7-9, 14]. 
Although significant progress in reducing the lattice errors is being made, many of the CKM matrix 
element determinations and most of the constraints on the CKM unitarity triangle are still limited 
by lattice uncertainties, and it will take several more years to obtain percent-level accuracy in many 
quantities. 

Current observations are consistent with Standard Model expectations that the CKM matrix is 
unitary to our present level of experimental and theoretical precision. Although we have not yet 
found a "smoking gun" of new physics in the flavor sector, there are several 2-3(7 hints that should 
be monitored, such as the tension between the £k band and the remaining CKM unitarity triangle 
constraints and the fy s puzzle. New physics may not appear first, however, in the gold-plated 
processes listed in Fig. 2 that are the focus of most lattice efforts. Thus, given the mature status of 
many lattice weak matrix element calculations, it is now time to expand the standard lattice QCD 
repertoire to include rare K- and 5-decays, beyond-the-Standard Model contributions to neutral 
kaon and 5-mixing, and other even more challenging quantities such as K —> kk decay, D°-D° 
mixing, and the neutron electric dipole moment. The lattice community must be prepared for new 
physics wherever it may arise in the flavor sector so that we can take advantage of this exciting 
opportunity to find physics beyond the Standard Model. 
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